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摘 要 : 塔 式 果 园 风 送 喷雾 机 是 目前 普遍 使 用 的 果园 喷雾 机 ， 塔 式 结构 易 产 生 旋转 和 不 规则 的 垂直 气流 ， 


导致 喷雾 机 出 风口 两 侧 风 场 分 布 不 对 称 ， 


分 布 规律 不 易 预 测 。 为 探究 适用 塔 式 果园 风 送 喷雾 机 出 风口 风 


场 建 模 方法 ， 本 研究 基于 喷雾 机 出 风口 风 场 计算 流体 动力 学 (Computational Fluid Dynamics, CFD) 建 模 方 


法 ， 提 出 了 用 


户 定义 函数 (User-Defined Function, UDF 


分 段 式 三 维 风速 入 口 边 界 条 件 设 置 方法 ， 并 研究 


了 淇 流 模型 和 计算 域 太 寸 对 喷雾 机 CFD 风 场 模拟 结果 的 影响 特性 。 采 用 Fluent 软 件 ， 建 立 了 三 种 CFD 模型 : 


模型 一 以 11 个 区 域 的 平均 风速 作为 边界 条 件 ; 模型 二 采 


步 研 究 计 算 域 尺寸 对 风 场 模拟 的 影 


J UDF 分 段 式 三 维 风速 入 口 作为 边界 条 件 ; 为 进 一 


自 ， 建 立 了 小 计算 域 尺寸 的 模型 三 。 三 种 模型 均 采用 基于 雷诺 时 均 控 制 


(Reynolds—Averaged Navier-Stokes, RANS) 方程 的 上 -s 消 流 模型 和 大 -wow 消 流 模型 进行 风 场 计算 。 为 了 验证 模 


型 的 可 靠 性 ， 设 计 了 空间 风 场 立体 测 


量 系 统 ， 实 现 了 精确 快速 空间 风速 测量 。 验 证 结果 表明 ，Standard k-e, 


Realizable k-e, BSL Kk-w 和 SST k-w@ 漠 流 模 型 较 适 合 风 场 CFD 建 模 ， 其 中 Standard k-e 满 流 模型 运算 结果 最 


优 ， 模 型 决定 系数 尼 为 0.81。 基 于 


研究 结果 可 为 喷雾 机 出 风口 风 场 CFD 建 模 方 法 提供 参考 。 
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用 格式 : 翟 长 远 , 张 燕 妮 , SED, 了 
(中 英文 ), 2021, 


边界 条 件 ; UDF; 济 流 模型 ;计算 域 


DF 分 段 式 三 维 风 速 入 口 边界 条 件 设置 方法 不 仅 提 高 了 仿真 结果 的 准确 
性 (提高 了 5%)， 而 且 降 低 了 计算 的 复杂 性 。 在 其 他 参数 设置 相同 的 情况 下 ， 大 尺寸 计算 域 模型 的 性 能 略 优 
于 小 尺寸 计算 域 。 实 际 建 模 过 程 中 建议 根据 计算 机 计算 


能 力 、 念 真 的 实际 要 求 设置 合适 的 计算 域 尺寸 。 本 
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1 引言 


果园 的 病虫害 防治 长 期 依靠 化 学 农药 ， 据 统 
计 ， 如 停止 使 用 化 学 农药 ， 果 园 将 损失 66% 一 
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匀 的 特性 °°, P ESR E T AE o 
B A Bi AR Be wh 25 FEE RAS, ADU BER 
药 、 污 染 环境 ， 更 对 果品 安全 生产 带 来 巨大 威 
胁 。 果 园 风 送 喷雾 风力 的 有 效 调控 是 降低 雾 滴 球 
移 、 提 高 农药 利用 率 的 关键 “ 。Landers ”通过 
试验 发 现 风力 调节 可 提高 30% 的 农药 沉积 、 减 少 
75% 药 液 漂移 。 风 力 调 控 主 要 是 通过 改变 喷雾 机 
风 箱 进出 风口 面积 、 风 箱 角度 、 风 机 转速 等 方法 
以 获取 最 佳 的 喷雾 风 场 “"。 

风 送 喷雾 机 在 不 同调 控 条 件 下 的 空间 风 场 模 
型 构建 是 进行 风力 精准 调控 的 前 提 ， 空 间 风 场 的 
获得 可 以 通过 田间 实际 测量 完成 。Holown- 
icki 等 1 对 关于 风 场 的 工作 进行 了 梳理 ， 研 究 者 
们 为 了 确定 风 场 参数 对 雾 滴 沉积 的 影响 ， 进 行 了 
大 量 的 风 场 试验 与 研究 ， 但 大 空间 多 风 场 测量 过 
程 时 间 和 经 济 成 本 巨大 ， 再 加 上 自然 风 方 向 和 大 
小 的 又 变 ， 田 间 试 验收 集 的 风 场 数据 ， 对 不 同 的 
环境 条 件 具 有 不 可 重复 性 。 近 年 来 ， 研 究 者 进行 
了 大 量 的 计算 流体 动力 学 (Computational Fluid 
Dynamics, CFD) 风 场 建 模 方法 研究 ， 模 拟 不 同 
FEF WE SS BLA Tb BUG, SER Ze CFD 风 场 模 
型 可 以 准确 模拟 不 同类 型 喷雾 机 的 外 部 风 
场 “7 。 对 于 特定 类 型 的 喷雾 机 ， 经 过 试验 数据 
验证 的 风 场 模型 将 是 一 种 便捷 、 可 有 效 分 析 喷 雾 
机 外 部 风 场 的 工具 “2 。 

CFD 风 场 建 模 主 要 包含 了 模型 初始 条 件 TE 
制 方程 与 汪 流 模型 的 选择 ， 以 及 设计 合理 的 计算 
域 3 个 方面 的 内 容 "。Salcedo 等 “采用 简化 的 
直线 代替 喷雾 机 出 风口 ， 采 用 基于 雷诺 时 均 控 制 
(Reynold-Averaged Navier-Stokes，RANS) 方程 
的 满 流 模型 建立 了 CFD 风 场 2D 模型 。Dekeyser 
等 引 采 用 和 矩形 截面 和 弯曲 截面 模拟 横流 式 风机 
和 轴 流 式 风 机 出 风口 ， 尽 可 能 靠近 出 风口 测量 出 
风口 风速 ， 采 用 基于 RANS 的 潮流 模型 建立 了 
CFD 风 场 3D 模型 。Endalenw 等 “测量 距离 出 
风口 15 cm 处 三 维 风速 作为 模型 的 速度 入 口 ， 采 
用 平面 和 曲面 组 合 的 几何 结构 模拟 不 同类 型 的 喷 
雾 机 出 风口 ， 采 用 基于 RANS 的 潮流 模型 构建 喷 


FA. Duga 等 “采用 和 矩 形 截面 模拟 出 风口 ， 
以 出 风口 风速 轮廓 作为 计算 的 初始 条 件 ， 使 用 非 
定常 雷诺 平均 方程 (Unsteady Reynold-Averaged 
Navier-Stokes, URANS) 和 Kk-e 模 型 求解 计算 域 
气流 。Garcia-Ramos 构 建 风机 的 内 部 结构 ， 以 进 
和信 风机 的 总 体 风 量 代替 出 风口 风速 作为 模型 的 人 
口 边 界 条 件 ， 使 用 ANSYS-CFX 软件 ， 基 于 
Spalart-Allmaras 消 流 模型 (用 于 求解 粘性 滑 流 运 
动 单方 程 模型 ) 建立 风 场 ”"。Hong 等 ” 以 圆柱 
体 代替 出 风口 ， 以 测量 出 风口 风速 值 作为 模型 的 
入 口 风速 值 ， 采 用 基于 RANS 的 SST find Ste E 
立 风 场 。Badules 等 ” 进行 出 风口 内 部 几何 结构 
建 模 ， 并 研究 了 不 同 汕 流 模型 下 仿真 结果 差异 ， 
表明 准确 的 出 风口 内 部 结构 建 模 显著 提高 了 仿真 
结果 的 准确 性 ， 选 择 合适 的 汕 流 模型 对 风 场 建 模 
至 关 重 要 。 

研究 者 针对 不 同类 型 喷雾 机 进行 了 出 风口 外 
部 CFD 风 场 建 模 方法 研究 ， 发 现 塔 式 果 园 风 送 
喷雾 机 将 出 风口 风力 与 药 量 由 垂直 方向 变 为 水 平 
方向 ， 可 提高 果树 冠 层 顶部 的 药 液 沉积 。 但 塔 式 
结构 内 部 导 流 装置 布置 不 合理 会 导致 出 风口 风速 
分 布 不 均匀 且 存 在 潮流 的 现象 。 为 探究 适用 塔 式 
果园 风 送 喷雾 机 出 风口 风 场 建 模 方法 ， 本 研究 采 
用 不 同 的 边界 条 件 和 汕 流 模型 ， 建 立 了 喷雾 机 空 
间 风 场 对 比 模型 ， 通 过 搭建 空间 风 场 立体 测量 试 
验 台 ， 进 行 了 风 场 测量 和 模型 验证 ， 以 获得 一 种 
不 同 边界 条 件 设置 方 法 ， 以 及 计算 域 尺寸 和 汕 流 
模型 组 合 的 最 佳 果园 喷雾 机 外 部 风 场 CFD 建 模 
方法 。 


2 CFD 风 场 建 模 


2.1 喷雾 机 结构 


试验 用 果园 风 送 喷雾 机 型 号 为 3WG-1200A， 
出 风口 最 大 高 度 204 cm。 喷 雾 机 风 箱 外 形 为 
塔 式 ， 风 箱 高 度 154 cm， 下 边沿 距 地 面 高 度 为 
50 cm， 出 风口 宽度 8.3 cm， 出 风口 处 安装 10 个 
导 流 板 ， 将 出 风口 分 为 11 个 格 段 ， 导 流 板 厚 度 
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Alem, QA 1 ta. 
2.2 控制 方程 与 建 模 方法 


喷雾 风 场 模型 采用 CED 商用 软件 Fluent16.0 


拖拉 机 药 箱 风 箱 喷头 出 风口 


图 1 试验 用 喷雾 机 
Fig.1 Sprayer for test 


建立 。 建 模 时 ， 假 设 空 气 是 恒温 、 不 可 压缩 的 粘 
性 牛顿 流体 ， 空 气 参数 为 经 验 值 (Fluent 的 默认 
值 )。 控 制 方程 包括 连续 方程 、 动 量 方程 和 注 流 
方程 ”。 密 度 恒定 的 定常 流动 的 连续 方程 如 公 
式 (1) 所 示 。 

du/dx + dv/dy + ðw/ðz = 0 (1) 

其 中 ,uu、v、w 是 速度 和 失 量 uw 在 x、y、z 方 问 
的 分 量 ，m/s。 

N-S 动 量 守 恒 方 程 用 于 描述 流体 运动 ， 但 其 
对 计算 机 内 存 和 性 能 要 求 非常 高 ， 本 人 研究 采用 基 
F RANS 动量 方程 的 庙 流 模型 模拟 喷雾 机 风 场 ， 
RANS 动量 方程 组 如 公式 (2) 所 示 。 


j =d; _ Op | alpu”) E a(pu'v’) E a(pu'w') 
div (puu) = div ( ugradu) ie F 7 | 
0 | ð uv’ ð aya 9 PE N. 
div (pvu) = div (ugradv) - > E uoo ) 人 ) ayy | 03 
i = di _ ap ,|_a(pu'w') _ a(pv'w') alow?) 
div (pwu) = div (ugradw) az | a a 二 
其 中 ,wu，v,，w 为 三 维 风速 的 平均 值 ，m/s; ”Prandtl 数 ;，S 和 5S, 是 用 户 定义 的 源 项 。 
oe de os pan ew) (5) 
E, kg/m’; p 为 气流 压力 ，Pa; S,. S,. S, 为 对 2 2 


MD, -puv Aw AI 〈 存 在 6 个 不 同 的 
雷诺 应 力 项 )。 

Kk-e 漠 流 模型 和 k-w 满 流 模型 是 风 场 建 模 中 
较 常用 的 油 流 模型 。k-e 湛 流 模型 的 输 运 方程 见 
公式 (3) 和 公式 (4), 

a(pk) . a(pku,) ð | a) a, | 

| “st + G, 4 


ot Ox; Ox; O) OX; 


G,=pe=Yy* 8, (3) 
d(pe) , A(peu;) ð fe e)a 


ot Ox; Ox; Cr) OX; 


E g 
Cult CG) Carts (4) 


其 中 ，4w 为 流体 的 动力 粘度 ，N s/m G, 
是 由 于 平均 速度 梯度 引起 的 济 动 能 的 产生 项 ; 
G, 由 于 浮力 引起 的 注 动 能 的 产生 项 ; Cus Ca 
和 Cs 为 经 验 和 常数 ; oo 和 5o, 分 别 是 与 注 动 能 ( 公 
式 (4) ) MERKE (公式 (6) ) 对 应 的 


其 中 ，w 为 时 均 速 度 的 张 量 符号 ，m/s; Kk 为 
mape, m/s? 

u, = pC, — (6) 

其 中 ，C,= 0.09， 为 经 验 常数 ; p 为 空气 密 
E, 默认 为 1.293 kg/m’; u, X m ai E, 
N - s/m’。 

对 于 定常 流动 ，0()/9t= 0。 不 可 压 流 体 由 
于 浮力 引起 的 G=0， 可 压 注 流 中 脉动 扩张 的 贡 
献 ,= 0， 且 在 方程 设置 中 用 户 无 自 定义 的 源 
项 ， 本 研究 中 上 和 s 的 输 运 方程 可 简化 为 : 


ee 2 k ea tG, -pe (7) 


Ox; Ox; Ox) Ox; 
0( peu,) 0 u,\ ðe E 
ð = ð KH ð | G,, k G, = 
Xj X Oy} OX; 
e 
Cap T (8) 
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KE W LD AY k-o Tia tat Be AY Yi HE k A EEFE 
BCR œ 的 输 运 方程 如 下 所 示 2 
a(pku;) ð Ee 


Ox; Ox; 


Jean (9) 


O, OX; 


d(pwu,) ð am 


Ox; OX; | Fy OX; 


| +P, 一 工 (10) 

w = elk (11) 

其 中 ，P, 和 P, 代 表 了 Kk 和 w 的 产生 项 ;7 和 
了 ,代表 了 涝 流 有 效 耗 散 项 ， 其 他 参数 与 清流 模型 
中 的 参数 定义 规则 一 致 。 

Standardk -es 模型 用 于 强 旋 流 或 带 有 弯曲 壁 
面 的 流动 时 会 出 现 失 真 。 RNG -se 可 以 更 好 地 处 
理 高 应 变 路 及 流 线 弯 曲 程 度 较 大 的 流动 。Realiz- 
able -es 模型 可 用 于 旋转 均匀 剪 切 流 ， 管 道内 流 
动 ， 边 界 层 流动 ， 以 及 带 有 分 离 的 流动 |。 
BLSk-w@ 济 流 模型 在 内 部 区 域 采 用 k-w id Tt 
型 ， 在 自由 剪 切 流 中 采用 丰 -s 庙 流 模 型 ， 避 免 了 
大 -中 模型 的 强 自由 流 敏 感性 。 考 虑 了 主 滑 流 剪 应 
力 的 传递 效应 ,在 BSLK-w 消 流 模 型 的 基础 上 
对 满 动 粘度 进行 修正 ， 即 SSTk-w 注 流 模型 3。 
为 探究 不 同市 流 模型 对 计算 结果 的 影响 ， 本 研究 
计算 了 Standard k-e, RNG k-e, Realizable k-e, 
SST k-@, BSL k-o Ail Standard k - œ 7x PP m i 
型 的 风 场 ， 并 进行 了 模型 对 比 。 

采用 半 隐 式 算法 (Semi-Implicit Method for 
Pressure-Linked Equations, SIMPLE) 进行 流 场 
数值 计算 ， 其 特点 是 在 交错 网 格 的 基础 上 求解 压 
力 场 ， 从 而 求解 动量 方程 ， 需 要 提供 初始 的 压力 
场 和 速度 场 ， 通 过 不 断 地 迭代 逐渐 逼近 压力 和 速 
度 的 真 值 。 离 散 格 式 为 二 阶 迎 风 ， 在 二 阶 迎 风格 
式 中 实际 上 只 是 对 流 项 采用 了 二 阶 迎 风格 式 ， 而 
扩散 项 仍 采用 中 心 差分 格式 。 模 型 中 ， 喷 雾 机 出 
风口 域 四 周 和 空 气 域 的 地 面 均 设 置 为 壁面 ， 壁 面 
粗糙 高 度 k = 0， 壁 面 粗 糙 度 常数 C, = 0.5, Tit 
流 强度 T= 5%, mi BELL w,/u = 10。 


2.3 边界 条 件 


为 对 比 不 同 边界 条 件 对 模拟 结果 的 影响 ， 


CFD 风 场 建 模 分 别 设置 两 种 边界 条 件 。 边 界 条 件 
一 设置 方法 : 喷雾 机 出 风口 被 10 个 导 流 板 分 制 
为 11 部 分 ， 对 被 导 流 板 分 割 的 出 风口 内 风速 平 
均 ， 生 成 图 2 (a) 所 示 散 点 ， 并 以 笛 卡 尔 坐标 输 
入 作为 模型 一 的 风速 入 口 。 

边界 条 件 二 如 图 2 (b) FR, 采集 的 风速 数 
据 使 用 C 语 言 写成 分 段 式 风速 方程 ， 再 采用 用 户 
定义 函数 (User-Defined Function, UDF) PA% 
FIFE KAR A oe Se A = SE a HZ, PE 
和 模型 三 风速 入 口 。 


ex =y 
200) et ay 2.00 = 
© z 
180 F 。 1.80 
Š 160 E is 
Ro b : im 
a 1.40 “: 2 1.40 
2 | : 2 gh 
i 120 上 。 à mo 120 
. 
il ae 1 
= 1.00 ‘ we = 1.00 
0.80 -。 id 0.80 
bed 
0.60 0.60 
0.40 0.40 


-2 -2 8 18 28 38 -12 -2 8 18 28 38 
测量 风速 / Cues!) WEKE (mest) 


(a)11 个 区 域 三 维 平均 风速 (b) 三 维 风速 曲线 
图 2 喷雾 机 出 风口 风速 


Fig. 2 Air velocity at outlet of the sprayer 


如 图 3 所 示 ， 计 算 域 分 为 空气 域 和 喷雾 机 出 
风口 ， 通 过 交界 面 连 接 空 气 域 和 喷雾 机 域 ， 风 力 
能 够 穿 过 交界 面 进入 空气 域 ， 四 周 设置 为 无 滑 移 
壁面 (图 3 (a) )。 空 气 域 (图 3 (b) ) AHRR 
去 交界 面 的 其 他 位 置 设置 为 压力 入 口 ， 压 力 为 0; 
侧面 及 顶 面 设 为 压力 出 口 ， 压 力也 为 0， 地 面 设 
置 为 无 滑 移 壁 面 。 


2.4 计算 域 与 网 格 划分 


如 图 4 所 示 ， 共 建立 三 种 对 比 模型 。 模 型 一 
(图 4 (a) ) 对 应 边界 条 件 一 ， 计算 域 10 m 
(W) X10m (H) X25m (L) (其 中 W 为 宽度 
width, 工 为 长 度 long，H 为 高 度 height); 模型 二 
(图 4 (b) )， 对 应 边界 条 件 二 ， 计 算 域 10 m 
(W) X10m (H) X25m (L)。 风 场 试 验 在 大 
田 进行 ， 即 开放 的 空间 。 为 了 研究 计算 域 斥 寸 对 
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Ath 


ghat 


10 m(H) 


a: aa) 
A 


JEJ tH ry 


边界 条 件 二 | be 
压力 入 口 交界 面 

注 : W 为 宽度 width, 工 为 长 度 long, HAH EK height 

图 3 喷雾 机 域 和 空气 域 条 件 设置 
Fig. 3 Computational domains consisting of atmosphere 
part and sprayer part 

仿真 结果 影响 特性 ， 建 立 了 模型 三 (4 (c) ), 对 

应 边界 条 件 二 , 计算 域 为 Sm(W)X7m (B) X 

10m (L)。 


y 


边界 条 件 


边界 条 件 - 模型 二 
(b) 模 型 二 示意 图 


WR || 模型 三 网 格 划分 细节 


(c) 模 型 三 示意 图 
图 4 关于 不 同 边 界 条 件 研 究 的 三 种 CFD 对 比 模型 
Fig. 4 Three comparison CFD models for different 


boundary conditions 


采用 ICEM CFD 16.0 建 立 几 何 结构 并 划分 网 
格 。 空 气 域 为 正六 面体 结构 网 格 ， 喷雾 机 域 为 四 
面体 非 结构 网 格 ， 网 格 数量 与 尺寸 参数 如 表 1 所 
示 。 模 型 一 由 于 其 几何 结构 ， 需 要 在 划分 网 格 的 


过 程 中 依据 不 同 的 出 风口 进行 分 块 划分 ， 除 此 之 
外 ， 三 种 模型 网 格 划 分 方法 相同 。 以 模型 二 为 
例 ， 在 喷雾 机 行进 方向 设置 230 个 网 格 ; 在 竖 直 
方向 ， 喷 雾 机 出 风口 网 格 尺 寸 为 0.04 m， 并 从 
0.04 m 逐渐 增 大 至 0.5 m; 在 下 风 方 向 ， 喷 雾 机 
出 风口 网 格 尺寸 为 0.04 m， 并 从 0.04 mm 逐渐 增 大 
到 0.5 m。 最 后 ， 细 化 喷雾 机 出 风口 附近 网 格 ， 
加 密 后 的 喷雾 机 出 风口 附近 网 格 几何 三 维 尺寸 为 
0.02m (W) x0.02m (H) x0.02m (L). 

#1 不 同 边界 条 件 研 究 三 种 模型 网 格 数量 与 尺 十 


Table 1 Grid parameters of three models 


网 格 数量 与 参数 设置 。 模型 一 模型 二 模型 三 
四 面体 网 格 数量 4201 4474 8242 
四 面体 网 格 参数 /m 0.02 0.02 0.01 
正六 面体 网 格 数 量 3,820,375 3,434,833 1,114,074 
正六 面体 网 格 参数 /m 0.02 0.02 0.02 


3 试验 材料 和 方法 


3.1 空间 风 场 立体 测量 试验 台 设 计 

针对 喷雾 机 空间 风 场 区 域 大 且 复 杂 难 以 快速 
采集 获得 的 问题 ， 搭 建 了 空间 风 场 立体 测量 试验 
a, WAS hia. 


Fe. 


att 制 器 


图 5 空间 风 场 立体 测量 试验 台 


Fig. 5 Three dimensional measurement test-bed for airflow field 


试验 台 由 移动 装置 步 进 电 机 控制 器 、 电 机 豫 
动 器 、 导 轨 、 滑 块 、 同 步 带 、 文 架 和 地 轮 、 风 速 
仪 安 装 支架 、 风 速 仪 ， 以 及 上 位 机 三 部 分 构成 。 
移动 装置 中 步 进 电机 控制 器 (KH-01 型 ,浙江 热 
龙 自动 化 有 限 公 司 ) 通过 驱动 器 (HB860D 型 ， 
深圳 市 鸿 宝 达 有 限 公 司 ) 精确 控制 步 进 电机 
(HB80D 型 ， 深 圳 市 鸿 宝 达 有 限 公司 ) 带动 滑 块 
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沿 导 轨 平 稳 移动 ， 滑 块 移动 速度 和 行程 精确 可 
控 。 风 速 仪 支架 与 滑 块 刚性 连接 ， 带 动 固定 在 其 
上 面 的 风速 仪 随 滑 块 精确 移动 。 

风速 仪 为 热 敏 式 风速 仪 (型 号 GM8903， 
深圳 市 聚 茂 源 科 技 有 限 公 司 ) ， 测 量 范 围 为 
0~30 m/s, IXH 0.001 m/s ， 采 样 频 率 为 
1 Hz， 与 上 位 机 通过 USB 连接 ， 并 用 Hot Wire 
Anemometer 软 件 实时 存储 并 显示 风速 值 。 


3.2 田间 试验 方法 


试验 在 北京 小 汤山 国家 精准 农业 示范 基地 开 
展 ， 喷 雾 机 拖拉 机 头 朝 向 正 西 ， 在 试验 过 程 中 测 
量 行进 方向 右 侧 的 风 场 。 

(1) 喷雾 机 出 风口 风速 测量 。 如 图 6 (a) 
所 示 ， 喷 雾 机 出 风口 按照 几何 形状 分 为 竖 直 面 、 
曲面 和 和 斜面 三 部 分 。 喷 雾 机 出 风口 测量 原理 和 过 
程 见 图 6 (b), 首先 沿 风 箱 外 边缘 用 马克 笔 以 
1 cm 测量 间距 标定 ， 风 速 仪 贴 着 标定 刻度 ， 在 
X、Jy、z 方 向 ， 于 喷雾 机 出 风口 左 、 中 、 碳 三 个 
位 置 ， 间隔 1 cm 测量 ， 共 测量 160 组 数据 。 将 同 
一 高 度 左 、 中 、 右 3 个 风速 值 取 平 均值 ， 得 到 
160 组 喷雾 机 出 风口 三 维 风速 。 


(a) 出 风口 测量 原理 图 (b) 试 验 现场 图 
图 6 喷雾 机 出 风口 风速 测量 

Fig. 6 Air velocity measuring test of outlet of sprayer 

(2) 空间 风 场 立体 测量 。 如 图 7 所 示 ， 在 喷 
Z TFAA 〈Y 方 向 ) ， 采 样 点 与 风 箱 出 风口 竖 
直 部 分 的 水 平 距离 (S) 分 别 为 03、0.5、1.0、 
15、2 和 3m (下 文 提 到 的 水 平 距 离 S 含 义 均 与 
此 处 相同 )。 在 竖 直 方向 〈Z 方 向 )， 取 10 个 采样 
高 度 ， 在 0.23 至 2.03 m 之 间 采 样 间隔 为 0.2 m。 
试验 过 程 中 ， 在 滑 块 去 程 进行 数据 测量 ， 滑 块 位 


移 为 3 m， 移 动 速度 为 0.02 m/s; 回程 不 进行 数 
据 测 量 ， 滑 块 移动 速度 为 0.3 m/s。 在 每 个 采样 高 
度 获 取 150 个 风速 值 ， 共 获取 60 组 数据 ， 可 以 覆 


ti PMR SS 28 TH] KU o 


DLE 


采样 起 始 位 置 SS ot 


注 :彩色 线条 代表 采样 位 置 
图 7 空间 风 场 测量 采样 分 布 图 
Fig. 7 Distribution of sampling locations of airflow field 
measurement 

喷雾 机 、 试 验 台 和 基准 台 相 对 位 置 如 图 8 所 
示 ， 激 光 尺 固定 于 基准 台中 心 。 首 先 ， 打开 激 光 
K (型 号 为 LM40， 优 利 德 科技 股份 有 限 公 司 )， 
调整 风速 仪 探 头 ， 使 激光 束 穿 过 探头 孔 中 心 CFL 
尺寸 为 0.6 cmX0.8 cm)， 确 保 每 一 次 测量 起 始点 
一 致 。 以 基准 面 水 平 度 为 基准 水 平 ， 使 用 水 平 仪 
(EM5416-150 型 ,乐清 市 伊 菜 科 电气 有 限 公 司 ) 
调整 试验 台 ， 确 保 试验 过 程 中 试验 台 的 俯仰 角 


图 8 喷雾 机 、 试 验 台 和 基准 台 相 对 位 置 图 
Fig. 8 Relative position map of sprayer, test-bed and 


base station 
超声 波 风速 风 回 传感器 安装 在 喷雾 机 侧 后 
方 ， 与 喷雾 机 前 进 方 向 夹 角 为 132"， 距 离 基准 面 
中 心 为 14.8 m， 安 装 高 度 2 m。 风 场 测量 试验 共 
持续 2 d， 试 验 过 程 中 自然 风 平 稳 ， 在 测量 过 程 
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中 自然 风 有 少量 时 间 超 过 3 ms， 但 整体 的 测量 
环境 符合 国际 上 风速 喷雾 机 作业 标准 ， 可 认 
为 田间 试验 风 场 测量 数据 可 靠 ， 作为 风 场 模 型 验 


模型 ， 使 用 计算 机 (Intel i5, RAM, 16GB, 
Win10) 计算 喷雾 风 场 。 
模型 网 格 数量 可 能 会 对 模型 精度 产生 影响 ， 


证 的 基础 数据 。 
4 仿真 结果 与 验证 


4.1 空间 风 场 立体 测量 试验 台 设 计 


如 表 2 所 示 ， 对 三 种 CFD 模型 采用 六 类 汕 流 


在 网 格 建 模 同时 进行 网 格 无 关 性 验证 。 对 于 模型 
一 和 模型 二 的 RNGK-s 模 型 ， 残 差 曲线 一 直上 下 
正弦 波动 ， 难 以 收敛 到 10 以下， 连续 性 残 差 偏 
大 。 模 型 二 Standardf-s 淇 流 模型 的 滑动 能 大 计算 
发 散 。 在 Badules 等 研究 中 ， 喷 雾 机 产生 气流 的 
复杂 性 使 得 连续 性 方程 的 残 差 难于 小 于 104， 当 


表 2 不 同 满 流 模型 计算 的 3 种 CFD 模 型 特性 


Table 2 Characteristics of three CFD models for different turbulence models 


te 模型 一 


模型 二 模型 三 


SST k-@,2.43 X 107 
BSL k-«@,2.47X 107 
Jina Ye Fo HEY Standard k-w,5.30 107 


SST k-@,2.70X 107 SST k-@ 4.05 X 107 
BSL k-@,9.71X 107 BSL k-@,2.31X10° 
Standard k- œ ,大 发 散 Standard k-¢,4.54X 107 


与 连续 性 残 差 Standard k-¢,4.09X 10° Standard k-¢, 1.91 107 Standard k-¢, 1.38 107 
RNG k-¢,3.30X 10"! RNG k-e, 1.8810"! RNG k-¢,9.12X 107 
Realizable k-¢,3.99 X 10° Realizable k-¢, 1.88 107 Realizable k-¢,2.76 107 
. RNG k-e iit Vit i ,连续 方程 残 Standard k- oo inh We EY, k FEB ,模型 不 收 eee 本 
运算 结果 ee ae, eee ee eh FA Bb HE DRE A, APAR 
aT He 差 偏 大 ,运算 结果 不 可 接受 ;其 余 GL RNG ke Hi a A , FELT PRR ZAMR 可 接受 
E z auy = <i EE mere Were y a 
m RAN FRI TY EA 于 107 ;其余 潮流 模型 运算 可 接受 
Velocity V 


30.51 


(a) 模 型 一 BSL k -w 


30.51 


(e) 模 型 二 BSLk-w 


30.51 


(i) 模 型 三 BSLk-w 


图 9 三 


(b) 模 型 一 SSTK-w 


(£) 4241 = SST k-« 


(j) 模 型 三 SSTk-@ 


(c) 模 型 一 Standard k-¢ (d) #4! — Realizable k -€ 


(g) 模 型 二 Standard k-e (h) 模 型 二 Realizable k -€ 


(k) 模 型 三 Standard k-e 
种 模型 中 心 面 风 场 分 布 云图 


(1) 模 型 三 Realizable k -€ 


Fig. 9 Distribution of airflow field on the center of three modes 
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EEN TEIRA cE TE 3x 107 IN, A A eae 
的 计算 结果 ， 计 算 中 止 H. ASE SST k - o, 
BSL k-@, Standard k-e ll Realizable k-e 四 种 沸 流 
模型 计算 风 场 的 连续 性 残 差 (运算 中 最 大 值 ) 均 
可 稳定 在 10*， 计算 结果 可 接受 。 

风 场 中 心 面 速度 最 大 ， 具 有 代表 性 ， 选 取 长 
为 Sm、 高 为 3 m 的 平面 ， 提 取 模 型 的 中 心 面 速 
度 云图 。 如 图 9 为 表 2 得 出 的 四 种 运算 可 接受 的 
满 流 模型 计算 的 喷雾 机 外 部 风 场 中 心 面 云图 。 模 
型 一 与 模型 二 ， 或 者 模型 一 与 模型 三 的 风 场 云图 
分 布 轮廓 存在 明显 差异 ， 可 见 边界 条 件 设置 方法 
是 模型 风 场 分 布 的 主要 影响 因素 。 对 比 模型 二 和 
模型 三 ， 可 以 发 现 计算 域 尺 对 风 场 有 微小 的 影 
响 。 对 于 同一 几何 模型 (以 模型 二 为 例 ) ， 当 汕 
流 模 型 不 同时 ， 风 场 分 布 也 有 微小 的 差异 。 

因此 ， 边 界 条 件 设置 方法 是 影响 风 场 分 布 的 
主要 因素 。 运 算 可 接受 的 满 流 模型 和 计算 域 的 尺 
才 是 影响 风 场 分 布 的 次 要 因素 。 


4.2 模型 验证 


图 10 为 水 平 距离 S=1.0 m， 距 地 高 度 H= 


1.03 m 的 实测 数据 曲线 ， 风 场 风 速 横向 分 布 呈 
现 不 完全 左右 对 称 的 山峰 状 ， 横 坐标 是 风 场 横 握 
距离 ， 风 速 仪 移动 速度 为 0.02 m/s， 采 样 总 长 度 


-0.3 -02 -0.1 0 01 02 03 
风 场 宽度 /m 


图 10 2% S=1 m,H=1.03 may KM) Maz wh 2 


Fig. 10 Measurement air velocity curve, at S=1 m,H=1.03 m 

75 SB SC ys UA FE EmA, ESN A ih 
线 峰 值 风 速 ， 如 表 3 所 示 ， 用 于 CFD 模 型 验证 ， 
当 H=1.03 m, S=2 mm 时 ,测量 位 置 风速 为 1.33 m/s, 
同一 高 度 处 S=3 m 时 ， 风 速 为 2.61 m/s。 一 般 来 
说 ， 距 离 出 风口 较 近 位 置 风速 应 大 于 较 远 位 置 风 
速 ， 可 能 由 于 自然 风 风 速 风向 变化 ， 使 得 实际 风 
场 分 布 与 理论 风 场 分 布 不 同 ， 类 似 的 点 还 有 了 H= 
0.63 m, S=1.5 m #ilH=2.03 m, S=1.5 m, 


表 3 喷 雳 风 场 不 同 高 度 和 不 同 水 平 距离 时 风 场 宽度 方向 的 实测 最 大 风速 


Table 3 Peak measurement air velocity of sampling points with different height and horizontal distances 


单位 :m/s 
水 平 距 离 /m 距 地 高 度 / m 

0.23 0.43 0.63 0.83 1.03 1.23 1.43 1.63 1.83 2.03 
0.3 1.099 5.962 14.73 22.03 17.65 21.07 11.78 15.7 11.36 7.208 
0.5 1.744 2.489 5.901 18.67 10.68 19.61 12.09 10.67 10.26 8.847 
1.0 1.151 5.362 4.866 8.034 11.13 10.5 12.72 9.014 6.39 8.544 
1.5 1.439 1.901 1.413 4.302 5.784 5.401 7.678 8.974 6.578 4.173 
2.0 2.126 1.604 1.932 2.554 1.33 3.336 5.845 6.593 2.717 6.337 
3.0 1.901 1.547 1.59 1.81 2.61 1.182 3.749 1.964 2.264 2.303 


在 模型 的 中 心 面 采 用 CFD Post 提 取 试 验 采 
样 点 位 置 的 风速 点 云 ， 结 合 表 3 数据 做 出 散 点 图 
和 趋势 线 图 ， 如 图 11 所 示 。 其 中 模型 二 对 应 
Standard k-e 漠 流 模 型 ， 实 测 风速 与 仿真 风速 散 
点 关于 拟 合 直线 聚合 性 最 佳 ， 模 型 的 决定 系数 
为 0.81。 空 间 风 场 试验 在 大 田 进 行 , 图 11 (g) 
标记 左下 角 的 红色 小 圈 内 ， 模 拟 风 速 值 均 小 于 仿 
真 风 速 值 ， 可 能 的 原因 是 在 喷雾 风 场 影响 微弱 的 


区 域 ， 上 自然 风 起 主要 作用 ， 使 得 该 区 域 测 量 风速 
大 于 模拟 风速 。 考 虑 到 田间 复杂 的 自然 环境 ， 风 
场 模型 的 决定 系数 R 仍 可 达到 0.81， 表 明 该 CFD 
风 场 模型 能 够 准确 预测 风 场 。 

模型 一 几何 建 模 过 程 复杂 ， 且 与 模型 二 相 
比 ， 其 模型 决定 系数 尼 俩 低 ， 说 明 采 用 一 个 完整 
的 出 风口 ， 通 过 UDF 导 入 风速 曲线 是 简单 准确 
的 建 模 方法 。 模 型 二 与 模型 三 相 比 较 ， 模 型 二 决 
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更 加 符合 实际 风 场 ， 但 计算 域 尺 寸 增加 ， 模 型 的 
运算 时 长 增加 ， 且 对 计算 机 性 能 要 求 也 大 幅 提 
高 。 大 尺寸 计算 域 的 优势 是 便于 后 期 进一步 模拟 
动态 的 、 整 个 果园 内 的 风 场 。 在 风 场 建 模 中 建议 
综合 考虑 运算 能 力 和 需要 建 模 的 果树 、 果 园 规 
模 ， 设 计 风 场 计算 域 。 
以 模型 二 为 例 (图 11 e) ~11 h) ), 不 
同 庙 流 模型 对 应 的 风 场 计算 结果 ， 其 散 点 分 布 与 
模型 决定 系数 存在 差异 。 表 明 在 CFD 风 场 建 模 
过 程 中 ， 选 择 适 合 的 潮流 模型 可 以 提高 模型 的 准 
确 性 。 
综 上 ， 以 喷雾 机 出 风口 三 维 风 速 分 段 折 线 作 
为 边界 条 件 ， 选 择 较 大 的 计算 域 和 Standard k - e 


场 最 接近 ， 是 最 佳 的 风 场 建 模 方法 。 
5 讨论 


由 喷雾 机 出 风口 风速 云图 (图 9) 可 见 ， 风 
速 变 化 复杂 ， 且 存在 满 流 。 如 图 2 (b), 喷雾 机 
出 风口 实测 风速 曲线 ， 在 不 到 20 cm 的 高 度 差 内 
风速 变化 超过 20 m/s。Z 向 风速 ， 在 大 于 1 m 
( 距 地 高 度 ) 时 ， 其 大 小 与 Y 向 风速 (垂直 出 风 
口 方向 ) 相近 ， 所 以 Z 向 风速 对 风 场 的 影响 将 不 
可 忽略 。X 向 风速 关于 轴线 左右 摆动 ， 这 是 造成 
实际 喷射 气流 左右 摆动 的 主要 原因 。 若 采用 一 维 
风速 作为 模型 的 风速 入口 ， 其 结果 将 不 能 反应 实 
际 中 复杂 湛 流 的 喷雾 机 ， 会 给 模型 带 来 严重 的 
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误差 。 

Hong 等 “使 用 4 个 风速 值 作为 速度 入 口 边 
界 条 件 却 得 到 了 和 较 好 的 仿真 结果 ,分 析 原 因 ， 其 
试验 用 喷雾 机 内 部 是 由 软 管 将 气流 分 离 ， 每 一 个 
出 风口 风速 较 稳 定 ， 而 本 研究 使 用 的 喷雾 机 出 风 
口 风速 分 布 极为 复杂 ， 且 存在 洋流 。 国 内 主流 风 
送 喷 雾 机 大 都 与 本 研究 喷雾 机 相似 ， 对 于 该 类 型 
的 喷雾 机 ， 无 法 照搬 Hong 等 的 建 模 方 法 ， 仅 考 
虑 出 风口 一 维 风 速 获得 的 仿真 结果 不 符合 实际 ， 
应 全 面 测 量 出 风口 三 维 风速 ， 以 相应 的 三 维 风速 
曲线 作为 模型 的 风速 入 口 边 界 条 件 。 

测量 喷雾 机 出 风口 三 维 风速 ， 借 助 UDF 生 
成 风速 曲线 ， 通 过 笛 卡 尔 坐标 导入 三 维 风速 ， 进 
ÍT CFD 风 场 模型 建立 ， 模 拟 结果 显示 ， 此 边界 
条 件 设置 方法 不 仅 降低 了 建 模 复杂 性 ， 而 且 具 有 
较 高 的 准确 性 与 工程 价值 。 边 界 条 件 是 决定 模型 
准确 性 的 关键 ， 采 用 UDF 出 风口 风速 折线 式 边 
界 条 件 设 置 方法 ， 输 入 三 维 风速 进行 风 场 建 模 是 
一 种 简单 准确 的 CFD 建 模 方 法 。 


6 4 it 


本 研究 通过 研究 风速 入 口 边 界 条 件 设 置 方 
法 、 计 算 域 和 汕 流 模 型 的 选取 对 CEFD 风 场 模型 
的 影响 ,确定 最 优 的 CFD 风 场 模型 建 模 方法 ， 
并 进行 了 田间 试验 验证 ， 结 论 如 下 : 

(1) 提出 了 果园 喷雾 机 风 场 CFD 模拟 UDF 
分 段 式 风速 入 口 边界 条 件 三 维 风 速 设置 方法 。 以 
Standard k-e m WREN H, KAH UDF IERA 
速 入 口 边界 条 件 (模型 二 ) 相 比 采用 离散 点 式 三 
维 风速 入 口 边 界 条 件 〈 模 型 一 ) ， 仿 真 风 场 准确 
性 提高 了 5%， 且 采用 UDF 分 段 式 风速 入 口 边界 
条 件 设置 方法 ， 几 何 建 模 简单 ， 模 型 运算 经 济 。 

(2) RM UDF 分 段 式 三 维 风速 入口 边界 条 
件 设置 方法 和 离散 点 式 三 维 风速 人 口 边 界 条 件 设 
置 方 法 , 10m (W) x10m (H) x25m (L) 
和 Sm (W) x7m (H) x10m (L)， 两 种 计 
算 域 尺寸 ， 建立 了 CFD 对 比 模型 。 验 证 试验 结 
RRHH, RH Standard k-e mM, ， 选 择 较 大 


的 计算 域 ， 仿 真 准 确 性 最 佳 ， 模 型 决定 系数 忆 
为 0.81。 

(3) 设计 了 风 送 喷雾 机 空间 风 场 立体 测量 系 
统 ， 以 2 cm 的 间隔 采集 风速 ， 获 取 10 个 高 度 、6 
个 平面 的 三 维 立 体 风 速 ， 具 有 定位 精确 、 测 量 快 
速 且 数据 完整 的 特征 ， 为 风 场 仿真 验证 试验 提供 
了 良好 的 测量 平台 。 


参考 文献 : 

[1] RAK, 杨 代 斌 , 同 晓 静 , 等 . 农药 有 效 利 用 率 与 喷雾 
技术 优化 四. 植物 保护 , 2011, 37(5): 14-20. 

YUAN H, YANG D, YAN X, et al. Pesticide effective 
utilization rate and spray technology optimization[J]. 
Plant Protection, 2011, 37(5): 14-20. 

何 雄 奎 . 中 国 精准 施 药 技术 和 装备 研究 现状 及 发 展 
ENDI. 智慧 农业 (中 英文 ) 2020, 2(1): 133-146. 
HE X. Research progress and developmental recom- 


[2 


a 


mendations on precision spraying technology and 
equipment in China[J]. Smart Agriculture, 2020, 2(1): 
133-146. 

FEAR wb, 赵 春江 , WANG N, 等 . 果园 风 送 喷雾 精准 控 
制 方法 研究 进展 四. 农业 工程 学 报 , 2018, 34(10): 
1-15. 

ZHAI C, ZHAO C, WANG N, et al. Research progress 
on precision control methods of air-assisted spraying in 
orchards[J]. Transactions of the CSAE, 2018, 34(10): 
1-15. 

MIRANDA-FUENTES A, RODRIGUEZ-LIZANA A, 
GIL E, et al. Influence of liquid-volume and airflow 


[3 


[n 


[4 


rae 


rates on spray application quality and homogeneity in 
super-intensive olive tree canopies[J]. Science of the 
Total Environment, 2015. 537: 250-259. 

LIAO J, LUO X, WANG P, et al. Analysis of the influ- 
ence of different parameters on droplet characteristics 


[5 


[n 


and droplet size classification categories for air induc- 
tion nozzle[J]. Agronomy, 2020, 10(2): ID 256. 
[6] YANG F, XUE X, CAI C, et al. Numerical simulation 
and analysis on spray drift movement of multirotor 
plant protection unmanned aerial vehicle[J]. Energies, 
2018, 11(9): ID 2399. 
LANDERS A J. Developments towards an automatic 
precision sprayer for fruit crop canopies[C]// Pitts- 


[7 


m 


burgh, Pennsylvania. Michigan, USA: American Soci- 
ety of Agricultural and Biological Engineers, 2010. 

李 龙 龙 , EE, 宋 坚 利 , 等 . 基于 变量 喷雾 的 果园 所 
动 仿 形 喷 筋 机 的 设计 与 试验 . 农业 工程 学 报 , 2017, 
33(1): 70-76. 

LI L, HE X, SONG J, et al. Design and experiment of 


[8 


[om 


80 


ChinaXiv 合 作 期 刊 


202302.00198v1 


chinaXiv 


[9] 


[10] 


[11] 


[12] 


[13] 


[14] 


[15] 


[16] 


[17] 


[18] 


[19] 


orchard automatic profiling sprayer based on variable 
spray[J]. Transactions of the CSAE, 2017, 33(1): 70-76. 
QIU W, ZHAO S, DING W, et al. Effects of fan speed 
on spray deposition and for targeting air-assisted spray- 
er in pear orchard[J]. International Journal of Agricul- 
tural and Biological Engineering, 2016, 9(4): 53-62. 
OSTERMAN A, GODESA T, HOCEVAR M, et al. Re- 
al-Time positioning algorithm for variable-geometry 
air-assisted orchard sprayer[J]. Computers and Elec- 
tronics in Agriculture, 2013, 98: 175-182. 

GU J, ZHU H, DING W. Unimpeded air velocity pro- 
files of an air-assisted five-port sprayer[J]. Transac- 
tions of the ASABE, 2012, 55(5): 1659-1666. 
SALCEDO R, PONS P, LLOP J, et al. Dynamic evalu- 
ation of airflow stream generated by a reverse system 
of an axial fan sprayer using 3D-ultrasonic anemome- 
ters. Effect of canopy structure[J]. Computers and Elec- 
tronics in Agriculture, 2019, 163: 1-14. 

陈 建 泽 , 宋 淑 然 , 孙 道 宗 , 等 . 远 射程 风 送 式 喷雾 机 气 
流 场 分 布 及 喷雾 特性 试验 . 农业 工程 学 报 , 2017, 33 
(24): 72-79. 

CHEN J, SONG S, SUN D, et al. Experiment on the 
distribution and spray characteristics of the air flow 


field of the long range air feed sprayer[J]. Transactions 
of the CSAE, 2017, 33 (24): 72-79 

HOLOWNICKI R, DORUCHOWSKI G, SWIECH- 
OWSKI W, et al. Variable air assistance system for or- 
chard sprayers; concept, design and preliminary test- 
ing[J]. Biosystems Engineering, 2017, 163: 134-149. 
DUGA A T, RUYSEN K, DEKEYSER D, et al. CFD 
based analysis of the effect of wind in orchard spray- 
ing[J]. Chemical Engineering Transactions, 2015, 44 
(2015): 289-294. 

HONG S W, ZHAO L, ZHU H. SAAS, a computer pro- 
gram for estimating pesticide spray efficiency and drift 
of air-assisted pesticide applications[J]. Computers and 
Electronics in Agriculture, 2018, 155: 58-68. 

BAHLOL H Y, CHANDEL A K, HOHEISEL G A, et 
al. The smart spray analytical system: Developing un- 
derstanding of output air-assist and spray patterns from 
orchard sprayers[J]. Crop Protection, 2020, 127: ID 
104977. 

ERIB, 祁 力 钧 , SSL HH Bos J PBL Pat SH HT Da Ma z 
Fe WL AA Ud BY CED 模拟 及 验证 四. 农业 工程 学 报 ， 
2015, 31(11): 46-53. 

WANG J, QI L, XIAQ. CFD simulation and validation of 
trajectory and deposition behavior of droplets around 


ry 


target affected by air flow field in greenhouse[J]. Trans- 
actions of the CSAE, 2015, 31(11): 46 -53. 

SALCEDO R, VALLET A, GRANELL R, et al. Euleri- 
an—Lagrangian model of the behaviour of droplets 


[20] 


[21] 


[22] 


[23] 


[24] 


[25] 


[26] 


[27] 


[28] 


[29] 


[30] 


produced by an air-assisted sprayer in a citrus or- 
chard[J]. Biosystems Engineering, 2017: 76-91. 

HONG S W, ZHAO L, ZHU H. CFD simulation of pes- 
ticide spray from air-assisted sprayers in an apple or- 
chard: Tree deposition and off-target losses[J]. Atmo- 
spheric Environment, 2018, 175: 109-119. 

ATAM E, HONG S W, ARTECONI A. Thermofluid 
modelling of large-scale orchards for optimal design 
and control of active frost prevention systems[J]. 
Energies, 2020, 13(2): 378. 

SALCEDO R, GRANELL R, PALAU G, et al. Design 
and validation of a 2D CFD model of the airflow pro- 
duced by an airblast sprayer during pesticide treat- 
ments of citrus[J]. Computers and Electronics in Agri- 
culture, 2015, 116: 150-161. 

DEKEYSER D, DUGA A T, VERBOVEN P, et al. As- 
sessment of orchard sprayers using laboratory experi- 
ments and computational fluid dynamics modeling[J]. 
Biosystems Engineering, 2013, 114: 157-169. 
ENDALEW A M, HERTOG M, GEBREHIWOT M G, 
et al. Modelling airflow within model plant canopies 
using an integrated approach[J]. Computers and Elec- 
tronics in Agriculture, 2009, 66(1): 9-24. 

ENDALEW A M, DEBAER C, RUTTEN N, et al. A 
new integrated CFD modelling approach towards air- 
assisted orchard spraying— Part II: Validation for dif- 
ferent sprayer types[J]. Computers and Electronics in 
Agriculture, 2010, 71(2): 137-147. 

DUGA A T, DELELE M A, RUYSEN K, et al. Devel- 
opment and validation of a 3D CFD model of drift and 
its application to air-assisted orchard sprayers[J]. Bio- 
systems Engineering, 2016, 154: 62-75. 
GARCIA-RAMOS F J, MALON H, AGUIRRE A, et 
al. Validation of a CFD model by using 3D sonic ane- 
mometers to analyze the air velocity generated by an 
air-assisted sprayer equipped with two axial fans[J]. 
Sensors, 2015, 15(2): 2399-2418. 

HONG S W, ZHAO L, ZHU H. CFD simulation of air- 
flow inside tree canopies discharged from air-assisted 
sprayers[J]. Computers and Electronics in Agriculture, 
2017, 149: 121-132. 

BADULES J, VIDAL M, BONE A, et al. Comparative 
study of CFD models of the air flow produced by an 
air-assisted sprayer adapted to the crop geometry[J]. 
Computers and Electronics in Agriculture, 2018, 149: 
166-174. 

王 福 军 .计算 流体 动力 学 分 析 , CFD 软件 原理 与 应 
用 [MI. 北京 : 清华 大 学 出 版 社 , 2004. 

WANG F. Computational fluid dynamics analysis, prin- 


ciple and application of CFD software[M]. Beijing: Ts- 
inghua University Press, 2004. 


ChinaXiv 合 作 期 刊 


202302.00198v1 


chinaXiv 


81 
[31] YANG Y, TSE K T, JIN X, et al. A numerical tree cano- and validation of a 2D CFD model of the airflow pro- 
py model and its application in computational wind en- duced by an airblast sprayer during pesticide treat- 
gineering simulation[C]// The 7th International Collo- ments of citrus[J]. Computers and Electronics in Agri- 
quium on Bluff Body Aerodynamics & Applications. culture, 2015, 116: 150-161. 
Amsterdam, Netherlands: Elsevier, 2012: 1429-1436. [34] BADULES J, VIDAL M , BONE A, et al. Comparative 
[32] MENTER F R. Two-equation eddy-viscosity turbu- study of CFD models of the air flow produced by an 
lence models for engineering applications[J]. AIAA air-assisted sprayer adapted to the crop geometry[J]. 
Journal, 1994: 1598-1605. Computers and Electronics in Agriculture, 2018, 149: 
[33] SALCEDO R, GRANELL R, PALAU G, et al. Design 166-174. 


CFD Modeling and Experiment of Airflow at the Air 
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Abstract: The tower-type sprayer produces swirling and irregular vertical airstream. The complex swirling results in airflow 
asymmetry between sides of the sprayer, and the vertical air velocity profile can be unpredictable when the rotational speed of 
the fan changes. The spray deposition is directly linked to the airflow pattern obtained from the sprayers. In order to study air- 
flow field of this type of air-assisted sprayer, a CFD (Computational Fluid Dynamics) model for the tower-type sprayer was de- 
veloped. A boundary condition setting method of UDF (User-Defined Function) sectional 3D air velocity was proposed. And the 
influences of turbulence models and the size of computational domain on CFD airflow simulation were studied. Using Fluent 
software, three different CFD models were established. The Model 1 took the average air velocity of 11 regions as the velocity 
inlet. The Model 2 used UDF segmented three-dimension air velocity line as the boundary condition. In order to further study 
the influence of the computational domain size on simulation, the Model 3 with a smaller computational domain was estab- 
lished. The turbulence model based on reynolds-averaged navier-stokes (RANS) control equation was used to calculate the air- 
flow field in all models. In order to verify the reliability of the model, a three-dimensional measurement system of airflow field 
was designed, which was used for accurate and fast velocity measurement. The results showed that the Standard k-e turbulence 
model, Realizable k-e turbulence model, BSL k-w turbulence model, SST k-w turbulence model were suitable, and the Standard 
k-e turbulence model was the best one. The CFD boundary condition setting method of UDF sectional three-dimension air veloc- 
ity could improve the accuracy of simulation, and reduce the calculation complexity. With the same settings of other parameters, 
the performance of the CFD model with larger scale calculation domain was slightly better than that with smaller computational 
domain. The size of computational domain should be set to the appropriate extent, considering the calculation capacity and prac- 
tical requirements of modelling. The research results could provide an important reference for CFD modeling of spray airflow 
field. 
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